Three Dimensional Ising Model, Percolation Theory and Conformal Invariance 
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The fractal structure and scaling properties of a 2d slice of the 3d Ising model is studied using 
Monte Carlo techniques. The percolation transition of geometric spin (GS) clusters is found to occur 
at the Curie point, reflecting the critical behavior of the 3d model. The fractal dimension and the 
winding angle statistics of the perimeter and external perimeter of the geometric spin clusters at the 
critical point suggest that, if conformally invariant in the scaling limit, they can be described by the 
theory of Schramm-Lowner evolution (SLE K ) with diffusivity of n = 5 and 16/5, respectively, putting 
them in the same universality class as the interfaces in 2d tricritical Ising model. It is also found 
that the Fortuin-Kasteleyn (FK) clusters associated with the cross sections undergo a nontrivial 
percolation transition, in the same universality class as the ordinary 2d critical percolation. 
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Conformal invariance plays the key role in the analyt- 
ical description of two-dimensional (2d) systems at their 
critical points [1]. One of the approaches to study the 
critical behavior of spin models such as g-state Potts 
model (q < 4) is the geometrical approach based on per- 
colation theory [2, 3] and also finite-size scaling (FSS) 
relations supplied by conformal field theory (CFT) [4, 5]. 
Fortuin-Kasteleyn (FK) bond representation of Potts 
model gives a close relation between percolation and ther- 
mal phase transition in which critical singularities can 
be represented in terms of FK clusters [6] . An FK clus- 
ter can be constructed from a "geometric" cluster — i.e., 
a set of nearest neighbor sites of like states, by assign- 
ing a bond between each pair of sites with a certain 
temperature-dependent probability. The FK clusters al- 
ways percolate right at the critical temperature T c . The 
percolation temperature for both FK and geometric clus- 
ters of a g-state Potts model coincide only in two dimen- 
sions. 

Critical behavior of the interfaces of these clusters in the 
scaling limit is conjectured, and in some cases proven, to 
be described by the theory of Schramm-Lowner evolution 
(SLE K ) [7]. The diffusivity k denotes the classification of 
conformally invariant interfaces whose value for the hull 
and the external perimeter (EP) of an FK cluster, k and 
k' respectively, satisfy a duality relation kk' = 16 [8] (for 
a review on SLE, see [9]). Another similar duality also 
holds in 2d, between the hull of an FK cluster and the 
hull of a spin cluster. Considering these two duality re- 
lations together, one finds that the hull of a spin cluster 
and the EP of an FK cluster can be described with the 
same k. Many connections between SLE and various 2d 
systems such as turbulence, spin glasses, growth models, 
sandpile models and disordered Potts model have been 
recently noticed and discussed [10]. 

For the Ising model, equivalent to q = 2 Potts model, 
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in three dimensions (3d) the percolation transition of ge- 
ometric spin (GS) clusters occurs at some temperature 
T p well below the Curie point T c [11]. Thus, contrary to 
the FK clusters, the 3d geometric clusters do not cap- 
ture the critical properties of the model. The question 
that arises is how the FK and GS clusters in a 2d cross 
sectional slice of the 3d Ising model behave and whether 
they reflect the critical properties of the 3d model. This 
is the main subject of the present letter. 

In this letter we present the results of extensive Monte 
Carlo study of 3d Ising model simulated by using Wolff's 
algorithm [12] based on single cluster update. We exam- 
ine the FSS hypothesis for some percolation observablcs 
on a 2d slice of the lattice, indicating that the percola- 
tion threshold of the geometric clusters coincides with 
the critical point. This is in agreement with the observa- 
tion of [13], in which the authors have found a power-law 
scaling behavior in the length distribution of the loops 
formed by an arbitrary cross section of the boundaries 
of 3d GS clusters at T = T c . In this paper, we under- 
take a detailed and systematic study of the statistical 
and percolation properties of all GS clusters on an arbi- 
trary cross section of the system. Note that disjoint GS 
clusters on a 2d cross section may belong to a 3d GS clus- 
ter. The scaling behavior of the GS clusters motivated us 
to examine possible conformal invariance of the critical 
geometric clusters which seem to undergo a continuous 
transition at the threshold. We find that the perimeter 
and the EP of a GS cluster at the critical point satisfy the 
duality relation and their fractal dimensions and winding 
angle statistics are compatible, in the scaling limit, with 
the family of conformally invariant curves i.e., SLE K with 
k = 5 and 16/5 respectively. This latter is in the same 
universality class as interfaces in tricritical Ising model 
in two dimensions. 

Our analysis of FK clusters on 2d slices also indicates 
that these clusters undergo a nontrivial percolation tran- 
sition at /3^ K /f3 c = 1.617(1), in the same universality 
class as the classical bond percolation. 

There exist some numerical investigations of 3d Ising 



2 



0.7 



0.6 



0.5 



0.4 




(P/P e -P/P e )L" 

0.998 0.999 1 1.001 1.002 

P/P c 




FIG. 1: (Color online) (a) The percolation probability P 3 of 
GS clusters, measured on 2d cross sections of 3d Ising model 
of different size L 3 , as a function of (3. Upper inset: FSS 
plots of the data with appropriate values v = 0.36 and p' c = 
1.001172,9c. Lower inset: the data collapse very close to the 
crossing point for two P s curves of sizes L = 40 and 70, with 
v — 0.36 and j3' c = 1.003316/3 C . (6) Size dependence of the 
crossing point /3' C (L, L') of different curves with the reference 
P s curve of L' = 40 (dotted-dashed line). 



model on special geometries with specific boundary con- 
ditions (be) [14, 15] which confirm a linear relation be- 
tween the scaling dimension of the operators of 3d sys- 
tems and their correlation lengths. This relation has been 
given by Cardy [16] by exploiting the conformal group 
transformations. 

In order to investigate the fractal structure of the GS 
clusters on 2d cross sectional slices of 3d Ising model, 
we first examine the FSS hypothesis for some percola- 
tion observables i.e., percolation probability P s , percola- 
tion strength M, and the mean cluster size \- F° r this 
purpose, the simulation of 3d Ising model is carried out 



by using the Wolff's Monte Carlo algorithm [12] on sim- 
ple cubic (SC) lattices of different size L 3 with free be, 
and sizes up to L = 150. After discarding the first L 3 
Monte Carlo sweeps for equilibration, we then analyze 
the spin configuration of a planar slice of the lattice par- 
allel to the xz — plane located at y = \ L/2\ ([--J denotes 
the integer part). We check whether it contains a GS 
cluster which spans the plane in just a certain direction 
e.g., z— direction. Therefore, we can obtain the proba- 
bility P s to have a spanning GS cluster at temperature 
T. We then measure the average mass of the largest 
spin-cluster, M. At each temperature, the averages are 
taken over 10 2 independent simulation runs and during 
each run another average is taken over 10 5 independent 
samples each of which is taken after L Monte Carlo steps 
[25]. This rather large number of averages is needed to 
obtain smooth curves close enough to the Curie point. 
The results are shown in Figs. 1 and 2 as a function of 
inverse temperature (3. 

As shown in Fig. 1, the curves P S {P) computed for differ- 
ent cross sectional lattice sizes for sufficiently near sizes, 
cross at the same point very close to the Curie point, 
implying that the scaling dimension of the percolation 
probability is zero. This agrees with what is expected 
from scaling theory [4] which states that the P s curves 
should have the form P 8 (fi) = P s [(/3/j3 c - \)L X I V \. The 
exponent v characterizes the divergence of the correla- 
tion length £ as the percolation threshold is approached, 
£ ~ \/3//3 c — Note that here £ is proportional to 

the spatial extent of the GS clusters in 2d slices of the 
3d Ising model [26]. 

By utilizing the data collapse, it is possible to esti- 
mate the values of the exponent v and the crossing point 
of the curves /3' c . To measure the quality of the collapse 
of the curves we define a function S(v,j3) as a function 
of the chosen values of v and j3 (the smaller S is indica- 
tive of a better quality of the collapse — see [19, 20] and 
appendix of [21]). We find its minimum S m i n ~ 1.24 for 
values v = 0.36(3) and (3' c = 1.0012(10)/3 C , confirming the 
percolation threshold of the 2d GS clusters at the Curie 
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FIG. 2: (Color online) (a) The average mass of the largest GS 
cluster M, and (b) the mean cluster size x as a function of /3. 
Inset: size dependence of \ a t P — Pc- 
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point. Applying these values for v and f3' c , rcscaling of 
the data for different lattice sizes fall onto a universal 
curve (see upper inset of Fig. 1(a)). However crossing 
point of the curves is slightly different from /3 C , but if 
one does not consider this discrepancy, the proper data 
collapse would not be achieved because, due to the rel- 
atively small value of v, the factor L x / U would magnify 
this discrepancy along the a;— axis and thus the curves 
do not collapse onto a single function. This effect will 
be more evident when one zooms in around zero on the 
x— axis. We find that the crossing point moves towards 
the Curie point for larger system sizes (Fig. 1(6)). Con- 
sidering the P s curve for L' = 40 as a reference curve, we 
find that the crossing point of the other curves for differ- 
ent sizes with this reference curve i.e., /3' C (L,L'), behaves 
like [(#(£, 40)/j8 c ) - 1] ~ L~ u , with uj = 1.02(5), which 
is independent of L'. Therefore, to obtain a proper col- 
lapse even close enough to the crossing point, we need to 
measure /3' c with high precision which is possible for only 
two P s curves. An example of such collapse is shown in 
the lower inset of Fig. 1(a), for system sizes L = 40 and 
70. However, our estimation for the exponent v is not 
accurate and has to be reevaluated by considering the 
corrections-to-scaling effects. 

In order to check that the coincidence between the per- 
colation threshold (3' c of the 2d GS clusters and the Curie 
point p c is not specific to the considered lattice geometry, 
i.e., SC structure, we repeated the Monte Carlo compu- 
tations for a face-centered cubic (FCC) lattice. Note that 
both f3' c and j3 c depend on lattice geometry. We find again 
their coincidence for the FCC lattice. 

We have also computed the percolation strength A/, 
as the order parameter in the percolation, and the mean 
cluster size x (Fig- 2). We find that, at the Curie point 
X scales as L 1 , with 7' = 1.46(4). These are indicative 
of a continuous transition of 2d GS clusters. 
These scaling arguments suggest a fractal structure for 
2d GS clusters, meaning that at (3 C there exists a set of 
GS clusters and their surrounding loops at all scales. 

In order to investigate manifestation of conformal in- 
variance in the critical behavior of the GS clusters on 
the cross sections, we analyze statistical properties of the 
boundaries of these clusters and compare them with the 
SLE curves. We first measure the fractal dimension of the 
cluster boundaries and check the duality relation. The 
fractal dimension of SLE curves is given by d/ = 1 + re/8, 
and the duality conjecture [8] states that the EP of SLE 
hulls for re > 4 looks locally as SLE curves but with a 
dual value k' = 16/ re. 

We undertook extensive Monte Carlo simulations with 
anisotropic geometry to obtain a reliable estimate of the 
fractal dimension at the Curie point j3 c (the value of j3 c 
has been determined to a high precision in Ref. [22]). We 
present the results of calculations for rather large system 
sizes of L x x L y x L z , with L y = L z = L, L x = 4L and 
10 < L < 290. The periodic boundary condition was ap- 
plied along all three directions. After equilibration, we 
collected a number of 10 5 (5 x 10 4 ) of independent spin 
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FIG. 3: (Color online) The average length of the perimeter 
(■) and EP (a) of a spanning GS cluster versus the system 
size L, compared with ones for SLE5 and SLE 16 / 5 curves (solid 
lines), respectively. The error bars are almost the same size 
as the symbols. 



configurations on a planar cross section of the 3d sam- 
ple located at y = [L/2\ for smallest (largest) system 
size. Each configuration was taken after the L x Monte 
Carlo steps. We then identified all GS clusters on each 
cross section of a fixed size L, and marked all spanning 
clusters (if exist) along the z— direction. For each span- 
ning cluster, a walker was applied to determine each of 
its perimeters which connect the lower boundary to the 
upper one along the z— direction. The walker was forced 
to start moving from a node on the dual lattice (which 
is also a square lattice) lied at the most left (right) site 
belonging to the spanning cluster on the lower bound- 
ary, by using the turn-right (left) tit-breaking rule [23] 
and thereby, for each size L, an ensemble of interfaces 
was obtained. It is shown in [23] that using these tie- 
breaking rules, one can obtain the interfaces on a square 
lattice without any ambiguity compatible with the con- 
formal invariant properties of critical lattice models. Ac- 
cording to each interface, the EP was also determined as 
the border of the corresponding interface after closing off 
all its boundary narrow passageways of lattice spacing a. 

For fractal interfaces, their average length I, measured 
in units of lattice spacing a (which is set to unity here) , 
scales as a function of the system size L as I ~ a(L/a) df 
where df is the fractal dimension. Applying this relation 
for the ensemble of the perimeters and the EPs, we find 
their fractal dimensions equal to df — 1.621(10) and dj = 
1.389(10), respectively. These results give rise to re = 
4.96(8) and re' = 3.11(8) for the perimeter and EP of a 
GS cluster, respectively, comparable with ones for SLE 
curves with diffusivity re = 5 and re' = 16/5, which satisfy 
the duality relation (we find rere' = 15.42(65), comparable 
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with the exact one rere' = 16). To show this consistency, 
we have compared the results of our simulations with 
those predicted for SLE curves in Fig. 3. The exponent 
5' = 1.23(1) obtained in [13] for the hull of the GS clusters 
is also in a good agreement with our result S = 2/df = 
2/(1 + re/8) ~ 1.23077 with re = 5. 

To verify further the consistency between the statisti- 
cal behavior of the interfaces of GS clusters in a 2d cross 
section of the 3d Ising model and the SLE curves with 
re = 5 and re' = 16/5, let us now examine their winding 
angle statistics which is predicted by the theory of SLE 
for critical interfaces [7]. We use the definition given in 
[24]. For each interface, an arbitrary winding angle is at- 
tributed to the first edge (which is taken to be zero here). 
The winding angle for the next edge is then defined as 
the sum of the winding angle of the present edge and the 
turning angle to the new edge measured in radians. It is 
shown that [24] the variance in the winding grows with 
the system size like 

(0 2 )=c+^lnZ, (1) 

where c is a constant whose value is irrelevant. As shown 
in Fig. 4, the winding angle statistics of GS interfaces are 
in a good agreement with those of the corresponding SLE 
curves. We find that re = 4.93(10) and re' = 3.15(10). 

The point which has to be explained and discussed 
here is the seeming discrepancy between the exponents 
obtained from the direct percolation properties of the GS 
clusters and the ones obtained by analyzing the fractal 
properties of their hull and external perimeter. The lat- 
ter study gives rise to a connection between statistical 
behavior of the boundaries of GS clusters in a 2d cross 
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FIG. 4: (Color online) Variance of the winding angle for the 
perimeter (■) and EP (a) of GS clusters on a 2d cross section 
of the 3d Ising model. The solid lines are set according to the 
Eq. (1) for comparison, with c ~ —1.37 and re = 5 for upper 
graph and with c ~ —0.84 and re = 16/5 for the lower graph. 
In the inset, the variance in semilogarithmic coordinates. 
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FIG. 5: (Color online) Percolation probability for FK clusters, 
Pf K , on a 2d cross section of 3d Ising model as a function 
of /3, for different lattice size L. Inset: the data collapse with 
v — 4/3 for 2d classical percolation, and f3^ K //3 C = 1.617. 



section of the 3d Ising model with those of FK clusters 
in a pure 2d tricritical Ising model, while the values of 
the geometrical exponents obtained in the first part of 
the paper are obviously not those predicted for tricriti- 
cal Ising model at its thermal critical point. Although 
we do not have a full description of such discrepancy, 
its resolution seems to lie on natural difference between 
GS and FK clusters. Moreover, the notion of conformal 
invariance in the statistical behavior of the hull and ex- 
ternal perimeter of GS clusters requires extraction of the 
statistics of the Loewner driving function and showing 
that it is a Brownian motion, which is very hard to do on 
small samples such as ours. These will be investigated in 
details in a future work. 

So far, we have focused only on the statistical prop- 
erties of GS clusters. We now present the results of the 
same study for the FK clusters built up on a 2d spin 
configuration by assigning a bond between each pair of 
nearest neighbor sites of the like spins with probability 
p = 1 — e~ 2 ^. Our analysis indicates that the percola- 
tion transition of the FK clusters occurs at f3^ K / j3 c = 
1.617(1) (shown by the vertical solid line in Fig. 5). 
This is obtained by computing the percolation proba- 
bility shown in Fig. 5, with the same procedure as for 
GS clusters. We find that the critical behavior of these 
clusters is in the same universality class as the 2d critical 
percolation. As shown in the inset of Fig. 5, all curves 
fall well onto a universal curve with the same correlation 
exponent i.e., v = 4/3 as for classical bond percolation. 
The percolation threshold of 2d bond percolation on a 
square lattice is located at p = 1/2 which is determined 
by the vertical dashed line in the Fig. 5, at (3 < , 
indicative of a nontrivial percolation transition at j3^ K . 
The latter is expected due to the condensation of the 
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majority spins at low temperatures. The analysis of the 
statistical properties of the perimeter and EP of FK clus- 
ters will appear in the forthcoming paper. 
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